Exact weak-coupling radius of the Holstein polaron in one, two, and three dimensions 



00 
OS 

o> 

< 



CO 

c3 



I 

C 

O 

o 



> ■ 

oo ■ 

-xf : 

m . 
oo ■ 
O ■ 

oo : 

On . 



c : 
o 

o ■ 



X 
S3 



Aldo H. Romero 1,3 , David W. Brown 2 and Katja Lindenberg 3 
1 Department of Physics, 
University of California, San Diego, La Jolla, CA 92093-0354 

Institute for Nonlinear Science, 
University of California, San Diego, La Jolla, CA 92093-0402 

3 Department of Chemistry and Biochemistry, 
University of California, San Diego, La Jolla, CA 92093-0340 
(February 1, 2008) 



We apply weak-coupling perturbation theory to the Hol- 
stein molecular crystal model in order to compute an electron- 
phonon correlation function characterizing the shape and size 
of the polaron lattice distortion in one, two, and three di- 
mensions. This correlation function is computed exactly to 
leading order in the electron-phonon coupling constant, per- 
mitting a complete description of correlations in any dimen- 
sion for both isotropic and arbitrarily anisotropic cases. Us- 
ing this exact result, the width of the polaron is character- 
ized along arbitrary directions. The width of the polaron 
thus determined disagrees in every dimension with some well- 
known characterizations of polarons, signalling in particular 
the breakdown of the adiabatic approximation and the char- 
acterizations of self-trapping associated with it. 

PACS numbers: 71.38. +i, 71.15.-m, 71.35.Aa, 72.90. +y 



As quantum quasiparticles, polarons (of whatever par- 
ticular variety) describe the states excitations correlated 
with the deformation or polarization quanta of a host 
medium. Often characterized in terms such as "large" 
and "small" and "free" and "self-trapped", perhaps the 
single most intuitively-accessible polaron property is its 
size. Our principal interest in this paper is the "po- 
laron radius" , which is, however, not a uniquely-definable 
quantity. In many approximate approaches, one encoun- 
ters polarons characterized by spatially-localized semi- 
classical wavefunctions for which a radius appears to be 
easily defined through the variance of the squared prob- 
ability amplitude of a localized electron. Such localized 
polarons are often described as "self-trapped". When 
such localized wavefunctions are not found, or are not 
found to be energetically favored, the polaron is often 
said to be delocalized or " free" and the polaron radius is 
then said to be infinite (i ll . 

These characterizations of polaron radii do not gener- 
alize to energy band theory where, owing to the trans- 
lational invariance of the crystal lattice, all states are 
delocalized. When a Bloch state is constructed from lo- 
calized functions, the " identity" of the localized function 
is lost and one can no longer characterize the properties 
of the quantum state in terms uniquely associated with 
the local function used to convey local correlation struc- 
ture into the delocalized quantum state [Q. These lo- 



cal properties are probed with correlation functions that 
reveal the internal structure implicit in the delocalized 
state |TJ|. 

In this paper, we focus on one such correlation func- 
tion that has been long and widely used to characterize 
polaron size in D dimensions: 
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This function can be viewed as measuring the shape 
of the polaron lattice distortion around the instanta- 
neous position of the electron, or, essentially equivalently 
in view of the strongly local character of the electron- 
phonon coupling, as an image of the electron density one 
associates with a localized polaron in semi-classical ap- 
proaches. Indeed, a common finding in the adiabatic 
approximation is that this correlation function can be 
expressed as the real-space autocorrelation function of a 
localized electron density, whose spatial variance a 2 is 
simply proportional to the square of the polaron radius 
R. In 1-D, one finds the particularly common result 
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that captures in one very simple relationship between the 
electron transfer integral J, the phonon frequency uj, and 
the electron-phonon coupling constant g much of what is 
commonly considered characteristic of polarons in 1-D. 

We assert that the 1-D characterization of polaron size 
as contained in (||) and the 2-D and 3-D characterization 
of polarons as distinguishably " free" or " self-trapped" fail 
to capture the size dependence of polarons in the weak 
coupling regime. We apply weak-coupling perturbation 
theory to derive the exact form of the correlation func- 
tion (P in 1-D, 2-D, and 3-D, and from this we derive the 
exact real-space variance that characterizes the polaron 
radius in any direction. This result is incompatible with 
adiabatic theory in 1-D, 2-D, and 3-D, but is compatible 
with results of other approaches not limited by the adi- 
abatic approximation. This incompatibility ultimately 
holds significant implications for the notions of polaron 
size and the concept of self-trapping generally. 

We use the Holstein Hamiltonian [jL4|Jl5| 



H — Hkin + H p h + Hi. 
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in which at creates a single electronic excitation in the 

rigid-lattice Wannier state at site n, and bt creates a 
quantum of vibrational energy in the Einstein oscillator 
at site ft. All sums are understood to run over the entire 
-D-dimensional lattice of edge-length N. Because there 
is no phonon dispersion in this model, and because the 
electron-phonon coupling is strictly local, H p h and H int 
are directly affected by the dimensionality and structure 
of the lattice only through the number of terms contained 
in each sum. It is in Hkin where lattice dimensionality 
and structure have their greatest influence; the Ji are the 
nearest-neighbor electronic transfer integrals along the 
primitive crystal axes, and the £i are unit vectors associ- 
ated with the primitive translations. We emphasize that 
the discrete, integer D-tuples n label primitive unit cells 
along the primitive axes a, 6, and c, such that the dis- 
crete vector analysis involved in this paper is formally in- 
dependent of lattice structure. As such, the above model 
encompasses all Bravais lattices, with the different lat- 
tice structures appearing only in the relative values of 
the three hopping integrals J a , and J c . 

When considering less than three dimensions, one may 
be legitimately concerned with possible distinctions be- 
tween a pure 2-D scenario, for example, and cases of 
"mixed" dimensionality wherein the electronic degrees 
of freedom may be essentially 2-D in layers of a 3-D 
medium. Within the scope of the results here presented 
there is no distinction between "pure" and "mixed" di- 
mensionalities so that the electronic dimensionality is the 
only relevant dimensionality. The 2-D case is therefore 
obtained by setting any one Ji equal to zero, and the 1-D 
case is obtained by setting any two Ji equal to zero. 

The core issue of concern in this paper is evident in 
a small but powerful way in the expected value of the 
interaction energy in the polaron ground state (k = 0), 
Ei n t = (0|-ffmt|0). This energy is related to the center 
value of the correlation function (|l|) through 
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in which e p is the small polaron binding energy in the ex- 
treme strong-coupling limit. Any non-zero value of Ei nt 
implies a non-zero value of Cq D ^ , and, since Ei n t is non- 
vanishing at any finite Ji and g, it necessarily follows 
that Cg D ' is non-zero in any interesting case. However, 
since the correlation function (nu satisfies the sum rule 



Y,?Cf ] = 1 for all Ji and g in any dimension, this non- 

trivially finite value of C D ^ implies the generic existence 
of non-trivial local character in the electron-phonon cor- 
relations that define polaron structure. This is concluded 
without appeal to any approximation method or spe- 
cific calculation; it is a general consequence of a non-zero 
electron-phonon interaction energy in the polaron model. 
This limited, but broad conclusion already appears to 
contradict some expectations for polaron structure in two 
and three dimension where infinite-radius "free" states 
are suggested to prevail in the weak-coupling regime. It 
is possible, in principle, for localized correlations to be 
compatible with an infinite polaron radius if the localiza- 
tion is sufficiently weak that its moments are all infinite. 
To fully settle this question, direct computation is re- 
quired. 

We now pursue Rayleigh-Schrodinger perturbation 
theory |ll|[l(|[l7j in the weak-coupling regime identify- 
ing the unperturbed Hamiltonian as Hq = Hkin + H p h 
and the perturbation as H' = Hi nt 
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General K-dependences can be obtained, but for our pur- 
poses we require only the k = results 
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where 



E(k) = -2^ J, cos h 
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is the free electron band energy at finite k. Consequently, 
one finds that 

Cf 1 = 2<^°> [D] | Cf 1 | [D]) + {g 2 } , (12) 

which to leading order in g (equivalent to order g 2 in the 
ground state energy) yields 
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which can be further developed into the appealing ex- 
pressions (with obvious generalization to D > 3) 
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In which I n (z) is the modified Bessel function of the first 
kind. Note that when any one Ji is set to zero, the con- 
tribution from the i axis is reduced to a delta function 
5 ri o, showing one way in which the 1-D and 2-D cases 
can be recovered straightforwardly from the 3-D case. It 
is also the case, however, that for any set of Ji, sum- 
ming C\t over one recovers C ? , and further summing 

the latter over either of the remaining ri recovers ; 
thus, the correlation function that results from averaging 
out the transverse dependences of a multi-dimensional 
correlation function results in the same 1-D correlation 
function that would have resulted from setting the trans- 
verse Ji equal to zero. Figure [l] shows several examples 
of the 1-D result (fb 
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FIG. 1. Weak-coupling 1-D correlation function C, 
Ja/huj = 1, 2, 4, and 8. 
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If instead of averaging out transverse dimensions one 
makes a 1-D cut through a selected portion of the corre- 
lation function (e.g. fb=f c = 0), one finds something 
other than 1-D character. Figure g illustrates the dimen- 
sionality dependence of the correlation function by com- 
paring 1-D transections of the correlation function in 1-D, 
2-D, and 3-D for an isotropic case. The 1-D cut through 
the D-dimensional correlation function is normalized by 
the value of that correlation function at r = 0, so that 
what is displayed is a fair relative measure of how rapidly 
correlations decay with radial distance from the center of 
the polaron. As might be anticipated from the form of 
( |i~4| ) - ([l6|), correlations are seen to decay more rapidly 
with increasing real-space dimensionality. 

Characterizing this multi-dimensional correlation func- 
tion in terms of a width measure involves selecting a di- 
rection along which to contract a variance tensor; thus, 
given an appropriate unit vector jl, one writes 
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FIG. 2. Correlation decay along the a axis in the 
weak-coupling limit in 1-D, 2-D, and 3-D as measured by 
C [ °\~JC { ° ] for J a = J b = Jc= Shu. 

The components of <r 2 can be computed from 

r 

where i, j — a, 6, or c, resulting in 
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where 
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Thus, along each of the primitive crystallographic axes, 
the real-space variance is simply proportional to the elec- 
tron transfer integral along that axis, and in a general 
direction is just the appropriate mixture determined by 
rotation. 

Contrary to much prevailing opinion, this result shows 
that in the weak-coupling regime: i) there are no signifi- 
cant qualitative or quantitative differences between 1-D, 
2-D, and 3-D polaron radii, ii) the polaron radius in 2- 
D and 3-D is not infinite, and iii) the polaron radius 
does not scale as J/g 2 hui as commonly expected, but as 
a/ J/Tilo. 

Since these are exact results of weak-coupling per- 
turbation theory, the only way for these conclusions 
to be invalidated would be for perturbation theory to 
fail dramatically at infinitesimal <?; however, there is no 
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evidence of such potential and catastrophic failure to 
be found. Weak-coupling perturbation theory is well- 
behaved in the vicinity of the Brillouin zone center where 
both the unperturbed and perturbed states are well- 
separated from the one-phonon continuum (|/t| < k c , 
where 2J max [l — cosk c ] = hoj). In 1-D, the results of 
weak-coupling perturbation theory are broadly validated 
by a variety of high-quality, non-perturbative techniques 
demonstrating the qualitative and quantitative accuracy 
of the former to non-trivially finite coupling strengths 
p8| p8[ . Results are more scarce in 2-D and 3-D; how- 
ever, where reliable non-perturbative results are available 
in the appropriate regime, they support rather than con- 
tradict the results of weak-coupling perturbation theory. 

On the other hand, methods invoking the adiabatic ap- 
proximation do not enjoy such extensive corroboration 
in the weak coupling regime. This has been explained 
through a demonstration that one of the key features of 
the adiabatic approximation - a locking relation between 
electron and phonon coordinates - eventually demateri- 
alizes as electron-phonon coupling is decreased at finite 

J/huj HJ. 

We are led to conclude from all these considerations 
that the adiabatic characterization of self-trapping and 
the self-trapping transition are strongly limited in 1-D, 
2-D, and 3-D, and break down definitively in the weak- 
coupling regime. 

This and related work suggests that the self-trapping 
phenomenon is qualitatively similar in 1-D, 2-D, and 3-D, 
and involves a smooth transition between distinct weak- 
and strong-coupling polaron structures having large and 
small radii, respectively. The strong-coupling structure 
is consistent with traditional notions of small polarons, 
including strong coupling perturbation theory and the 
adiabatic approximation. The structure of polaron states 
in the weak-to-intermediate coupling regime appears to 
differ significantly from what has traditionally been asso- 
ciated with the terms "large polaron" and "free" , but still 
appears to be beyond the reach of many traditional meth- 
ods. A number of relatively new and mutually-consistent 
non-perturbative approaches to polaron theory afford re- 
liable and accurate information in these areas. We have 
used the Global-Local variational method to analyze 1-D 
polaron structure in detail throughout the polaron pa- 
rameter space p5L58j, and higher-dimensional calcula- 
tions by such methods are in progress pCfl . 
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